Chiral Perturbation Theory and UnitarizationQ 

E. Ruiz Arriola*, A. Gomez Nicola**, J. Nieves* and J. R. Pelaez** 

* Departamento de Ffsica Moderna, Universidad de Granada, E-18071 Granada, Spain. 
**Departamento de Ffsica Teorica Universidad Complutense. 28040 Madrid, Spain. 

Abstract 

We review our recent work on unitarization and chiral perturbation theory both in the 
7T7T and the irN sectors. We pay particular attention to the Bethe-Salpeter and Inverse 
Amplitude unitarization methods and their recent applications to irir and irN scattering. 
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O ■ 1 Introduction 

Chiral Perturbation Theory (ChPT) is a practical and widely accepted effective field theory 
to deal with low energy processes in hadronic physics. |l], g]. The essential point stressed 
in this approach is that the low energy physics does not depend on the details of the short 
distance dynamics, but rather on some bulk properties effectively encoded in the low energy 
\q , parameters. This point of view has been implicitly adopted in practice in everyday quantum 

physics; well separated energy and distance scales can be studied independently of each other. 
The effective field theory approach also makes such a natural idea into a workable and systematic 
computational scheme. 
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The great advantage of ChPT is that the expansion parameter can be clearly identified a 
priori (see e.g. Ref. [||) when carrying out systematic calculations of mass splittings, form factors 
and scattering amplitudes. However, the connection to the underlying QCD dynamics becomes 
obscure since the problem is naturally formulated in terms of the relevant hadronic low energy 
degrees of freedom with no explicit reference to the fundamental quarks and gluons. In addition, 
in some cases (see below) a possible drawback is the lack of numerical convergence of such an 
expansion when confronted to experimental data, a problem that gets worse as the energy of the 
process increases. Recent analysis provide good examples of both rapid convergence and slow 
convergence in ChPT. In the tttt sector the situation to two loops @ seems to be very good for 
the scattering lengths. Here, the expansion parameter is m 2 / (An f^) 2 = 0.01 (m v = 139.6MeV 
the physical mass of the charged pion) and the coefficients of the expansion are of order unity. 
For instance, to two loops (third order in the expansion parameter) the expansion of the s-wave 
of the isospin 1 = channel reads 

aoomTr = 0.156 + 0.043 + 0.003 + 0.015 + 0.003 + • • • (1) 

tree i l oop 2 loops 

where the theoretical errors are described in |J. Thus the expansion up to two loops is both 
convergent >> a^ n+l ^ and predictive Aa^ « a(™ +1 ) , Aa^ « s™. The prediction of 
ChPT of s-wave scattering lengths au in the isospin 1 = and 1 = 2 channels yields || 

aoom^ = +0.214 + 0.005 (exp. + 0.26 ± 0.05) 

a 20 m 7T = -0.420 + 0.010 (exp. -0.28 + 0.12) (2) 



1 Talk given by E.R.A. at the Miniworkshop on "Few Quark Problems", Bled, Slovenia, July-2000. 
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The theoretical predictions for these observables are an order of magnitude more accurate than 
the corresponding experimental numbers. Note that in hadronic physics we are usually dealing 
with the opposite situation, confronting accurate measurements to inaccurate theoretical model 
calculations. On the contrary, for irN scattering, for ChPT in the heavy baryon formulation 
the expansion is less rapidly converging than in the irir case, since NLO corrections become 
comparable to the LO ones. For instance in the i^-channel the expansion up to third order, 
after a fit to the threshold properties, reads 

a\ z ml= 353 + 4^95 - L49 +••• = 81.8 ± 0.9 (exp.80.3 ± 0.6) (3) 

1st order 2nd order 3rd order 

Here first order means 1//^, second order 1/ f^M^ and third order 1//^ and Despite 
these caveats, there is no doubt that the effective field theory approach provides a general 
framework where one can either verify or falsify, not only bulk properties of the underlying 
dynamics, but also the dynamics of all models sharing the same general symmetries of QCD. 

Finally, let us remark that the perturbative nature of the chiral expansion makes the gener- 
ation of pole singularities, either bound states or resonances, impossible from the very begining 
unless they are already present at lowest order in the expansion. There are no bound states 
in the irir and irN systems, but the p and the A resonances are outstanding features of these 
reactions dominating the corresponding cross sections at the CM. energies y/s = m p = 770MeV 
and \/s = 771a = 1232MeV respectively. 



2 The role played by unitarity 

Exact Unitarity plays a crucial role in the description of resonances. However, ChPT only 
satisfies unitarity perturbatively. Nevertheless, there are many ways to restore exact unitarity 
out of perturbative information, i.e.: the if-matrix method Ref. 0, the Inverse Amplitude 
Method (IAM) Ref. §, l, the Bethe-Salpeter Equation (BSE) Ref. pi, 111, the N/D method 



Ref. |iOjl , etc. ( See e.g. Ref. ]13[ for a recent review), which are closely related to one another. 



Some of these Unitarization methods have been very successful describing experimental data 
in the intermediate energy region including the resonant behavior. Despite this success the 
main drawback is that this approach is not as systematic as standard ChPT, for instance in 
the estimation of the order of the neglected corrections. In this work we report on our most 
recent works related to unitarization, both in the irir and ttN sectors, where we have obtained 
intermediate energy predictions using the chiral parameters and their error bars obtained from 
standard ChPT applied at low energies. This means in practice transporting the possible cor- 
relations among the fitting parameters obtained from a fit where the phase shifts are assumed 
to be gauss-distributed in an uncorrelated way. 



3 Results in the tttt sector 



The ChPT expansion displays a very good convergence in the meson-meson sector. As a conse- 
quence the unitarization methods have been very successful in extending the ChPT applicability 
to higher energies. In particular, within a coupled channel IAM formalism, it has been possible 
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Figure 1: Diagrammatic representation of the BSE equation. It is also sketched the used kinematics. 



to describe all the meson-meson scattering data, even the resonant behavior, below 1.2 GeV 
§, but without including explicitly any resonance field. These works have already been exten- 
sively described in the literature and here we will concentrate in the most recent works of two 
of the authors on the tttt sector, dealing with the Bethe-Salpeter equation, which has the nice 
advantage of allowing us to identify the diagrams which are resumed. 

Indeed, any unitarization method performs, in some way or another, an infinite sum of 
perturbative contributions. At first sight, this may seem arbitrary, but some constraints have to 
be imposed on the unitarization method to comply with the spirit of the perturbative expansion 
we want to enforce. In the BSE approach the natural objects to be expanded are the potential 



and the propagators. The BSE as it has been used in Refs. 11, 12] reads (See Fig. [y 



lj,(p,k) = V£(p,k) + if -^T P (q,k)A(q + )A(q^V^ P ,q) 



(27T) 



(4) 



where q± = (P/2±q) and Tp(p, k) and Vp(p, k) are the total scattering amplitude^ and potential 
for the channel with total isospin I = 0, 1, 2. and then the projection over each partial wave J 
in the CM frame, Tjj(s), is given by 



Tu(s) = ~ Pj (cos 6) T P (p, k) d(cos 9) 



i87TS 



Aa (s, m 2 ,m 2 ) 



1 



(5) 



where 9 is the angle between p and k in the CM frame, Pj the Legendre polynomials and 
A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. Notice that in our normalization the unitarity limit 
implies \T u (s)\ < 16tts/A 1 / 2 (s , m 2 ,m 2 ). 

The solution of the BSE at lowest order, i.e., taking the free propagators for the mesons and 
the potential as the tree level amplitude can be obtained after algebraic manipulations described 
in detail in Refs. [|ll|, [l^j , renormalization and matching to the Taylor expansion up to second 
order in s — Am 2 of the one loop chiral perturbation theory result. We only quote here the result 
for the p channel: 
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(7) 



2 The normalization of the amplitude T is determined by its relation with the differential cross section in the 
CM system of the two identical mesons and it is given by da/dQ. = \T P (p, fc)| 2 /647r 2 s, where s = P 2 . The phase 



of the amplitude T is such that the optical theorem reads ImTp(p,p) = — otot 

(s 2 - 4sm 2 ) 1/2 , with 

cross section. The contribution to the amputated Feynman diagram is (— iTp(p, k)) in Fig. |l|. 



the total 



3 



Figure 2: I — J = 1 irn phase shifts as a function of the total CM energy y/s for both sets of Z's given in 
Eq. (^). Left (right) figures have been obtained with the set A (B) of parameters. Solid lines are the predictions 
of the off-shell BSE approach, at lowest order, for the different I. J— channels. Dashed lines are the 68% confidence 
limits. Circles stand for the experimental analysis of Refs. and MM. 



where the unitarity integral Io(s) reads 

I (s) = Io(a) - / (4m 2 ) = -J—Ji-—l g \ 8 — (X) 




Here the complex phase of the argument of the log is taken in the interval [— 7r,7r[. Similar 
expressions hold for the scalar-isoscalar (a) and the scalar-isotensor channels. Notice that in 
this, so-called off-shell scheme, the left hand cut is replaced at lowest order by a pole in the 
region s « 0. For the low energy coefficients £1,2,3,4 we take the values 

set A : h = -0.62 ± 0.94, l 2 = 6.28 ± 0.48, l 3 = 2.9 ± 2.4, l A = 4.4 ± 0.3 

set B : h = -1.7 ± 1.0 , l 2 = 6.1 ± 0.5 , l 3 = 2.9 ± 2.4, l 4 = 4.4 ± 0.3 (9) 

In both sets Is and I4 have been determined from the SU(3) mass formulae and the scalar radius 



as suggested in [[jj and in 14], respectively. On the other hand the values of l\^ come from 
the analysis of Ref. |l5[ of the data on K^— decays (set A) and from the combined study of 
if/4— decays and tttt with some unitarization procedure (set B) performed in Ref. [O]. The 
results for the p channel are shown in Fig. y. 



As discussed in Ref. [12] it is also possible to take into account the left hand cut in the so- 
called on-shell scheme where it can be shown that after renormalization the on-shell unitarized 
amplitude acquires the following form 

Tjjis)- 1 +Io(s)-V I j(sr 1 = Tjjisor 1 +Io(s )-V IJ (s )- 1 = -C u (10) 

where Cu should be a constant, independent of s and the subtraction point so, and chosen to 
have a well defined limit when m — » and 1// — ► 0. Vu(s) is the on-s/ieZ/-potential and has the 
important property of being real for < s < 16m 2 , and presenting cuts in the four pion threshold 
and the left hand cut caused by the unitarity cuts in the t and u channels. The potential can 
be determined by matching the amplitude to the ChPT amplitude in a perturbative expansion. 
This method provides a way of generating a unitarized amplitude directly in terms of the low 
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energy coefficients £1,2,3,4 and their errors. In this on-shell scheme a successful description of 
both 7T7T scattering data as well as the electromagnetic pion form factor, in agreement with 
Watson's theorem, becomes possible yielding a very accurate determination of some low energy 
parameters. The procedure to do this becomes a bit involved and we refer to Ref. [12] for further 



details. We should also say that the one-loop unitarized amplitudes generate the complete ChPT 
result and some of the two and higher loop results. The comparison of the generated two-loop 
contribution of threshold parameters with those obtained from the full two loop calculation is 
quantitatively satisfactory within uncertainties. 



4 Results in the irN sector 

The methods and results found in the tttt system are very encouraging, suggesting the extension 
to the ttN system. However, ChPT does not work in the ttN sector as nicely as it does in the 
7T7T sector. As we will see, the low convergence rate of the chiral expansion makes it difficult to 
match standard amplitudes to unitarized ones in a numerically sensible manner. After an initial 
attempt within the relativistic formulation [19], it was proposed to treat the baryon as a heavy 



particle well below the nucleon production threshold [20|. The resulting Heavy Baryon Chiral 
Perturbation Theory (HBChPT) provides a consistent framework for the one nucleon sector, 
particularly in ttN scattering ||. The proposal of Ref. [^] adopting the original relativistic 
formalism but with a clever renormalization scheme seems rather promising but unfortunately 
the phenomenological applications to irN scattering have not been worked out yet. 



4.1 The I AM method in 7r-N scattering 

The inverse amplitude method (IAM) is a unitarization method where the inverse amplitude, 
and not the amplitude, is expanded, i.e., if we have the perturbative expansion for the partial 
wave amplitude f(u>), 

/M = / 1 M + / 2 H + / 3 M + ... (11) 

with co = y/q 2 + m 2 and q the CM. momentum, then one considers the expansion 

1 1 /2M+/3M , /2M 2 | _ (12) 



The IAM fulfills exact unitarity, Im/(a;) _1 = — q and reproduces the perturbative expansion to 
the desired order. 



This method has been applied for irN scattering [22| to unitarize the HBChPT results of 
Ref. § to third order. In this context, it is worth pointing out that the use of a similar unita- 
rization method, together with very simple phenomenological models, was already successfully 



undertaken in the 70's (see [23] and references therein). Nonetheless, a systematic applica- 
tion within an effective Lagrangian approach was not carried out. In ]^2| the phase shifts for 
the partial waves up to the inelastic thresholds have been fitted, obtaining the right pole for 
the A(1232) in the P33 channel. In that work, it has been pointed out that to get the best 
accuracy with data, one needs chiral parameters of unnatural size, very different from those 
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of perturbative HBChPT. This is most likely related to the slow convergence rate of the ex- 
pansion. However, it must be stressed that one can still reproduce the A(1232) with second 
order parameters compatible with the hypothesis of resonance saturation . In a subsequent 
work [^4|] we have proposed an improved IAM method based on a reordering of the HBChPT 
series. The encouraging results for the A-channel have also been extended to the remaining low 
partial waves [25], as it can be seen in Fig. § In this case, the size of the chiral parameters is 
natural and the x 2 P er d.o.f is considerably better than the IAM applied to plain HBChPT. 
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Figure 3: Phase shifts in the IAM method as a function of the CM. energy yfs. The shaded area corresponds to 
the result of propagating the errors of the chiral parameters obtained from low-energy data (see Ref. P5|). This 
illustrates the uncertainties due to the choice of different parameter sets from the literature. The dotted line is 
the extrapolated HBChPT result. The continuous line is an unconstrained IAM fis to the data, whereas for the 
dashed line the fit has been constrained to the resonance saturation hypothesis. 



4.2 BSE method and the A-resonance 



Recently f2q| , we have used the BSE to HBChPT at lowest order in the chiral expansion and have 
looked at the P33 channel. We have found a dispersive solution which needs four subtraction 
constants, 

ff/lx^r 1 = - { J^ m2) + PH + (^ 2 - ^VoH/ej 

PH = m 3 f Co + Cl (^-l)+ C2 (^-l) 2 + C3 (--l) 3 ) (13) 
V m m m ) 



where the unitarity integral is given by 



Jo(w) = Jo(w) — Jo(m) = — yL °. r, m {arcosh— — iir}; uj > m (14) 

4tt z m 
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Figure 4: P33 phase shifts as a function of the total CM energy y/s. The upper solid line represents a x 2_n t of 
the parameters, Co, 1,2, 3 to the data of Ref. [^lj (circles). Best fit parameters are denoted Co* lj2 ,3 in the main text. 
The lower lines stand for the results obtained with the parameters deduced from HBChPT and denoted Cg 1 2 3- 
Central values lead to the solid line, whereas the errors on 0^1,2,3 lead to the dash-dotted lines. 

The x 2 fit yields the following numerical values for the parameters: 

cf = 0.045 ± 0.021 , cf = 0.29 ±0.08, cf = -0.17 ± 0.09 , cf = 0.16 ± 0.03 

with x 2 /d.o.f . = 0.2. However, if we match the coefficients with those stemming from HBChPT 
we would get instead the following numerical values: 

cf = 0.001 ± 0.003 , cf = 0.038 ± 0.006 , cf = 0.064 ± 0.005 , cf = 0.036 ± 0.002 

The discrepancy is, again, attributed to the low convergence rate of the expansion. The results 
for the P33 phase shift both for the fit and the MonteCarlo propagated errors of the HBChPT 
matched amplitudes have been depicted in Fig. |] 

4.3 tvN scattering and the iV*(1535) resonance 

One of the greatest advantages of both the IAM and the BSE methods is that the generalization 
to include coupled channels is rather straightforward. For the IAM case we refer to || for more 
details. As for the BSE, two of the authors have dealt with the problem in the Sn channel, 
at y/s up to 1800GeV using the full relativistic, rather than the heavy baryon formulation used 
in Ref. [^J for the s-wave and extended in Ref. |28| to account for p-wave effects. For these 
energies there are four open channels, namely ttN ,r]N ,KT, and KA, so that the BSE becomes a 
4x4 matrix equation for the I = 1/2, J = 1/2 and L = partial wave. Additional complications 
arise due to the Dirac spinor structure of the nucleon, but the BSE can be analytically solved 
after using the above mentioned off-shell renormalization scheme. It turns out |3(J that to lowest 
order in the potential and the propagators, one needs 12 unknown parameters, which should be 
used to fit experimental data. Several features make the fitting procedure a bit cumbersome. 
In the first place, there is no conventional analysis in the relativistic version of ChPT for this 
process and thus no clear constraints can be imposed on the unknown parameters. Secondly, the 
channel 7riV — * irnN is not included in our calculation. Therefore one should not expect perfect 
agreement with experiment, particularly in the elastic channel since it is known that 10 — 20% of 
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Figure 5: nN scattering BSE results as a function of CM. energy \fs. Upper left figure: Sn phase shifts. 
Upper right figure: inelasticity in the nN channel. Lower left figure: nN — » r/N cross section. Lower right figure: 
7riV — > KA cross section. Data from Ref. Q. (See Ref. ||] for further details.) 



the N* resonance decay width goes into irnN. On the other hand, one cannot deduce from here 
how important is the ttttN channel in the r\ production channel, ttN — ► 7]N. Actually, in Ref. [29] 
it has been suggested that the bulk of the process may be explained without appealing to the 
three body intermediate state ttttN. The work of Ref. [29] would correspond in our nomenclature 
to the on-shell lowest order BSE approximation, which by our own experience describes well the 
bulk of the data. With the BSE we have provided a further improvement at low energies by 
including higher order corrections. In the absence of a canonical low energy analysis it seems 
wiser to proceed using the off-shell renormalization scheme. In Fig. || we present a possible 12 
parameter fit which accounts both for the elastic low and intermediate energy region and the 
lowest production channels. The failure to describe data around the N* resonance is expected, 
since as we have already mentioned the ttttN channel must be included. Our results seem to 
confirm the assumption made in Ref. [29] regarding the unimportance of the three body channel 
in describing the coupling of the N* resonance to rjN. 



5 Conclusions and Outlook 



The results presented here show the success and provide further support for unitarization meth- 
ods complemented with standard chiral perturbation theory, particularly in the case when res- 
onances are present. But unitarization by itself is not a guarantee of success; the unitarization 
method has to be carefully chosen so that it provides a systematic convergent and predictive 
expansion, as we have discussed above. In order to describe the data in the intermediate energy 
region the chiral parameters can then be obtained from 
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• Either a direct x 2 fit of the order by order unitarized amplitude and the corresponding 
low energy parameters. The upper energy limit is determined by imposing an acceptable 
description x 2 /DOF ~ 1. 

• Or from a low energy determination of the low energy parameters with errors by performing 
a x 2 -fit of the standard ChPT amplitude until x 2 /DOF ~ 1, and subsequent MonteCarlo 
error propagation of the unitarized amplitude. 

Differences in the low energy parameters within several methods should be compatible within 
errors, as long as the Chiral expansion has a good convergence. But, unfortunately this is not 
always the case. Clearly, the irN sector is not only more cumbersome theoretically than the 
7T7T sector but also more troublesome from a numerical point of view. Standard ChPT to a 
given order can be seen as a particular choice which sets higher order terms to zero in order 
to comply with exact crossing but breaking exact unitarity. The unitarization of a the ChPT 
amplitude is also another choice of higher order terms designed to reproduce exact unitarity but 
breaking exact crossing symmetry. Given our inability to write a closed analytic expresion for an 
amplitude in a chiral expansion which simultaneously fulfills both exact crossing and unitarity 
we have preferred exact unitarity. This is justified a posteriori by the successful description of 
data in the intermediate energy region, which indeed suggests a larger convergence radius of the 
chiral expansion. 
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